Optimum Placement of Post-IPN GW Chirp Templates 
Made Simple at any Match Level 
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A simple recipe is given for constructing a maximally sparse regular lattice of spin-free post-lPN 
gravitational wave chirp templates under a given minimal-match constraint, using Tanaka-Tagoshi 
coordinates. Cardinal interpolation among the resulting maximally-sparse templates is discussed. 
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I. INTRODUCTION 



H 

Gravitational wave chirps emitted by compact binary stars in their adiabatic inspiral phase will be primary targets 
for the early operation of broadband laser- interferometric detectors like TAMA300, GEO600, LIGO and VIRGO [Q. 
These signals have been thoroughly studied and accurately modeled 

The maximum-likelihood strategy J3| for detecting signals of known shape (except for a set of unknown parameters) 



in additive stationary colored gaussian noise consists in correlating the data with a set C = {gi} of possible expected 
^v^j \ waveforms (templates), and using the largest correlator as a detection statistic ||. The mentioned correlators are the 
. (noise-weighted) scalar products: 



(1.1) 



where (fi,f s ) is the antenna spectral window, a(f) — h(f)+n(f) are the noise-corrupted (spectral) data, h(j) is a 
(possibly null) signal, n(f) is a realization of the antenna noise with (one-sided) power spectral density n(/), g(f) is 
a template, and c.c. denotes complex-conjugation. 

The issue of optimum placement of the templates in the waveform parameter-space is more or less obviously a 
H ' crucial and still open one, and has been addressed by several Authors Q-Q- 

. At fixed false-alarm probability, the detection probability is an increasing function of the overlap 



between the signal h and the template g, where the norm ||u|| = (ujit) 1 / 2 . The template set C should be designed 
in such a way that, for any admissible signal h, the overlap never drops below a prescribed value T, which can be 
immediately related to the fraction (1 — T 3 ) of potentially observable sources which might be missed out ||. 

The overlap can be readily maximized w.r.t. the (unknown but irrelevant) initial-phase and coalescency-time 
differences between the signal and the template Q . The resulting partially-maximized overlap is called the match, is 
denoted by M(h,g), and is a function of the source and template intrinsic parameters only (e.g., at 2.5PN order, the 
binary companion-masses, spin-spin and spin-orbit parameters). The minimum-overlap condition rephrases obviously 
into the following minimal-match prescription: 

WieS, BgeC : M(h,g)>T, (1.3) 

where S is the set of admissible signals. 

Discussion will be henceforth restricted to the 2D parameter-space of spin-free chirp- waveforms, in view of the 
widespread present consensus that these should be adequate for the early operation of interferometric detectors pi . 

In the following we shall denote as G the (2D, spin-free) parameter-space point corresponding to waveform g. We 
shall also denote as jr(G) the match contour-line, whose points represent the set of waveforms h such that M(h, g) = T, 
and as Sr{G) the 2D-region bounded by jr(G), representing the set of waveforms h such that M(h,g) > T. 

The template set A = {Gi} should fulfill the following basic requirements. First of all, (a) the minimal-match 
condition ([hq) should be met, which re-phrases into 

(j5 r (G l )^S, (1.4) 
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where E is the image of S i.e., the subset of the waveform parameter-space corresponding to admissible sources. 
Pictorially, (|lJ) means that the patches S'r(Gi) should cover the whole set E without leaving holes. At the same 
time, (b) the template set A should be chosen as much sparse as possible, so as to minimize the overall number of 
templates, and hence the detection threshold-level, for any prescribed false-alarm probability ||. 

It might be further desirable (c) to deal with a set A forming a regular (or piecewise-regular) grid, for computational 
ease. 

A further template-placement requirement is also worth mentioning. When covering E with (regularly spaced 
or patch-wise regularly spaced) templates, a certain amount of template spill-over beyond 9E is unavoidable. The 
post-lPN spin-free waveform parameter-space subset E corresponding to admissible sources, for which m min <mi < 
rn 2 < m max , is a three-vertex curved-side 2D domain. Spill-over across the mi =m 2 boundary line of E (equal-mass 
boundary-line) is a true penalty, and thus (d) should be kept to a minimum. This is not the case for the other two 
boundary lines of E, where points beyond <9E could well correspond to possible (though unlikely) sources. 

At large values of the minimal match, e.g. typically M(h,g) > T > 0.97, as prescribed for a single-step search 
strategy ||, the match is very well approximated by a quadratic function JjJ of the distance between H and G in 
parameter-space, and the (approximate) match contour-lines 7r(G) are ellipses. However when using post-lPN order 
waveforms (which will be required in a one-step search to achieve the prescribed large minimal- match levels ||), as 
G moves throughout E, the above ellipses rotate and stretch as an effect of the intrinsic nonzero curvature of E. 

An effective procedure for constructing a 2D spin-free post-lPN template set appropriate to this case has been 
formulated in M , as a generalization of the 1PN strategy introduced in 0] . The simplest rectangular-cell template- 
lattice is chosen, and the span of template Gi (i.e., the set of points representing waveforms for which M(h,gi) > 
M(h, gj),Vj 7^ i) is accordingly taken as the largest rectangle inscribed in the minimal-match contour-line 7r(Gj). 
In order to minimize spill-over across the equal-mass boundary- line of E, while optimizing template placement for 
equal-mass sources (which are credited to be most abundant JToj] ) , the algorithm is started by placing a sufficient 
number of templates along the equal-mass boundary line. The nearest-neighboring templates are placed in such a 
way that the boundaries of their rectangular spans touch without intersecting. As a result E is covered with a set of 
patches, each of which is a stack of equal- width rectangles, whose height changes according to the shear/stretch of 
the elliptical 7r contour-lines. The current release of GRASP |Tl| ] implements a similar procedure, resulting into a 
relatively straightforward template placement algorithm. 

As already noted in |8) the above algorithm is not optimal, i.e., it does not necessarily satisfy requirement (b), in 
view of its a-priori restriction to rectangular cells []l2[ . 

The template placement issue becomes even more complicated at low minimal-match levels (e.g., typically r>0.7) 
as prescribed in the first step(s) of hierarchical search strategies. In multi-step hierarchical searches a single template- 
lattice A designed to achieve a large T throughout E is constructed. In the first step only the correlators corresponding 
to a suitably decimated subset of A are computed, covering E at a lower minimal-match level. In the second step, only 
a limited number of correlators, corresponding to a few non-decimated subsets of A need to be actually computed, 
covering suitable neighbourhood(s) of the candidate signals discovered in the first step. Multi-step hierarchical search 
strategies should allow a sizeable reduction in the required computing power |l3|-|l^]. Mohanty |l4| formulated a 
simple strategy for evaluating the lattice decimation-steps (taken as constant throughout extended patches of the 
lattice) for a post-lPN template lattice. Also in this case, however, requirement (b) will be most certainly violated. 
By itself, the template placement issue at low minimal-match levels is further complicated by the fact that the match 
contour-lines 7r(G) besides varying with G are no longer elliptical. For this reason, the procedure expounded in |T^| 
can only be validated by extensive Monte-Carlo trials. 

In this rapid communication we present a possible way to get rid of all the above difficulties and limitations, 
by capitalizing on the remarkable properties of the post-lPN waveform parametrization introduced by Tanaka and 
Tagoshi 16 , and further discussed in p7fl . We accordingly introduce a simple and systematic procedure to set up an 
optimum template-lattice Jig) ], which closely fulfills all conditions a), b) and c) above, at any prescribed minimal- 
match level [ jl9| . 

The paper is organized as follows. In Sect. 2 the Tanaka- Tagoshi parametrization and its relevant properties, 
including the main features of the match contour-lines, are briefly introduced. In Sect. 3 our optimum template 
placement procedure is exposed in detail. In Sect. 4 the gain achievable by cardinal interpolation among these 
optimally- placed templates is discussed. Conclusions follow under Section 5. 
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II. TANAKA-TAGOSHI PARAMETRIZATION OF POST-1PN WAVEFORMS 

The Tanaka-Tagoshi transformation p(| is a linear one-to-one mapping F 6 £ <-» F 6 T between the (2D, spin-free, 
post-lPN) admissible waveform parameter-space £ and a globally flat manifold T such that: i) the ratio 

_ M(h,g)-M(h,g) 

M{h,g) (2A > 

is negligibly small (typically of the order of 10~ 3 ]T^|) and, ii) M(h,g) depends only on the distance between H and 
G Jl6| . The Tanaka-Tagoshi construction is effective for all first generation interferometers [[Tt) . 

Clearly, property i) above allows to solve the optimum template placement problem in the globally flat manifold T, 
by finding a set A = {G,} C T such that: 

\JS r (G t )DT. (2.2) 

i 

On the other hand, property ii) implies that the sets Sr{Gi) are identical up to trivial translations. This means that 
A C T is a globally uniform template lattice, featuring a regular grid of nodes, at any prescribed minimal-match level. 

A. Constant Match Contour-Lines 
in Tanaka-Tagoshi Coordinates 

The match contour-lines 7r(G) in the (dimensionless) Tanaka-Tagoshi (spin-free) waveform parameter-space 
coordinates (xi,X2) have different geometrical features, depending on the value of T. In Fig. 1 the inverse of the 
curvature radius p 7 of 7r(G) is displayed as a function of the polar angle <fi around G for several representative values 

of r, for LIGO-I at 2.5PN order. The contour-lines 7r are also shown in the insets, where Sxt = x\ — x\ . 

In the limit as T — > 1 the contour-lines 7r are asymptotically circular and p~ x is almost constant. See, e.g., the 
r = 0.99 contour-line in Fig. 1 top-left. 

In a range r* < T < 1, the contour-lines 7r deviate from circular, but still remain globally convex. For LIGO-I, 
r* w 0.9225. Correspondingly, p^ 1 remains non negative throughout a full rotation around G. See, e.g., the T = 0.93 
contour-line in Fig. 1 top-right. 

For r < r*, the 7r contour-lines are no longer globally convex. In a range T' < T < T*, they consist simply of two 
convex [^lj arcs and two concave ones joining smoothly, to form a "peanut" shape. See e.g. the T = 0.914 contour-line 
in Fig. 1 bottom-left. For LIGO-I, w 0.9109. As T is decreased below T', the contour-lines become more and more 
complicated, as more and more bumps and dents do appear. Correspondingly, the boundary curvature p~ x exhibits 
more and more turning points throughout a full rotation around G. See, e.g., the T = 0.875 contour-line in Fig. 1 
bottom- right. 

For all values of T, the surface Sr(G) has always a centre of symmetry in G. 

III. OPTIMUM TEMPLATE PLACEMENT 

In the following we shall discuss the optimum template placement strategy for each of the above mentioned 
minimal-match ranges. Optimum here and henceforth means that requirements (a) , (b) and (c) discussed in Sect. 
1 will be strictly fulfilled JT|. 

A. The Asymptotic r — ► 1 Limit 

In the asymptotic limit r — > 1, the match contour-lines are circular. Rotational symmetry of the match contour- 
lines implies that the templates should sit at the vertices of regular polygons in the waveform parameter-space. These 
regular (open) polygons should besides make up a (regular) tiling of the (Tanaka-Tagoshi) plane f22|| , and hence can 
only be triangles, squares or hexagons. Let U be one such polygon, and {Gj} its vertices. The sparsest templates 
satisfying ( |2.2| ) are such that the circles 7r (Gi) touch at a single point P, which is obviously the center of symmetry 
of U. This means that the curve 7r(-P) goes through all the GVs. Hence, the templates should be located at the 
vertices of a regular triangle, square or hexagon inscribed in a circle 7r- 
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The obvious question is which polygon is best in terms of the total number of templates needed to cover the 
admissible-waveform parameter-space T in the Tanaka-Tagoshi plane. Clearly, the best polygon will be the one for 
which the span of each (and any) template has the largest measure (area) [p3| . 

The span (or, in more technical language, the Voronoi set j24j) of Gi is the set of points H for which M(h,gi) > 
M(h,gj), Vj ^ i, which are the images in the Tanaka-Tagoshi plane of those waveforms h which will be detected 
using template gi. The Voronoi sets are shown shaded in Fig. 2, for the (regular) triangular, square and hexagonal 
tilings of the (Tanaka-Tagoshi) plane pjj . 

A choice can be made on the basis of a comparison among the measures (areas) of the pertinent Voronoi sets. Let 
//(•) denote the measure (area), the p-gonal tiling-cell, and the corresponding Voronoi set. It is readily seen 
that M: 



/./, 



1.299037, 



y(6)j 



0.649519, 



(3.1) 



The triangular tiling is thus seen to be the best choice, yielding a gain in terms of template-span of ~ 30%, whereas 
the hexagonal one features a loss of ~ 35%, w.r.t. the simples square tiling p7| . 

The triangular tiling yields a template lattice whose fundamental domain (lattice-cell) is a parallelogram jl8| . Only 
three vertices of any lattice-cell belong to a single curve 7r ■ 

B. The Range T > V* 

In the range T > T*, the contour-lines 7r(G) are no longer circles, but still globally- convex a nd c enter-symmetric. 
Circular symmetry being lost, we shall still seek the sparsest template collocation subject to ( |2.2| ) among the non- 
regular, but still center-symmetric (open) p-gons (p = 3,4,6) with all vertices on yr, which tile the Tanaka-Tagoshi 
plane. 

It makes sense to compare for this case the template-span (Voronoi set) measures (areas) n[V (p) ], p = 3, 4, 6 to the 
template-span measure (area) ^[^o^l of the square tiling |2(| pertinent to the naive (but usual) approximation of 
circular match contour-lines. The ratios 



(4)i 



p = 3,4,6. 



(3.2) 



are displayed in Fig. 3 as functions of T for LIGO-I. Not unexpectedly, the triangular tiling is once more the best. 
This conclusion is valid for all first generation antennas, and the pertinent values of have been collected in Table-I. 



Antenna 


f3 


TAMA300 


1.58 


LIGO-I 


1.43 


GEO600 


1.58 


VIRGO 


1.38 



Table I.- The ratio r$, eq. (|3.2|), for first generation interferometers at T = 0.97. 



C. The Range T < T* 



By an argument of continuity (i.e., since the match contour-lines corresponding to T < T* can be obtained from 
those corresponding to T > T* by continuous transformations) we shall still seek the optimum tiling of Tanaka-Tagoshi 
plane using the (open) triangles with (all) vertices on 7p. However, the match con tour -lines being no longer globally 
convex, the largest-area triangle inscribed in 7r does not necessarily cope with eq. (2^), as e.g. exemplified in Fig. 4. 

One should accordingly seek thre e points Gi on 7r yielding the largest-area triangle G0G1G2, subject to eq. (|2.2|). 
The minimal-match condition (2.2) is equivalent to 



C 



|J [cnS r (Gi) 



(3.3) 



5 



wher e G is the (parallelogram) fundamental-domain of the lattice corresponding to a given triangular tiling. Condition 
(3.3) is most easily checked, involving only four patches S-p- 

Let P, P' ,Q,Q' the points (ordered, e.g., counterclockwise) where 7r and its convex-hull [^8| detach, shown in Fig. 
5. One can readily prove that two vertices of the sought optimum triangular-tile, say Go, Gi, should be sought on 
the arc PP' (or QQ') of 7r, while the third should be sought on the arc QQ' (resp. PP'). It is also readily proved 
that the lattice vector GqGi, after a trivial translation making its midpoint coincident with G, should be contained 
in the butterfly-shaped region (PGQ' U QGP'^j n S r (G) depicted in Fig. 5. 

The maximum-area triangular tiles subject to the constraint ( |2.2| ), are shown in Fig. 6, for several values of T in the 
range 0.7 < T < 1 (note that different scales are used in the the figure, for better readability). The lattice-cell measure 
corresponding to the optimum triangular-tile is displayed in Fig. 7 as a function of T, together with the number of 
templates needed to cover the companion-mass range O.2M < mi < mi < IOMq, for a number of representative 
minimal-match values. 



IV. CORRELATOR BANK ECONOMIZATION VIA CARDINAL INTERPOLATION 

In recent papers it has been shown that a sensible reduction (ss 75% for 1PN and higher-PN-order templates) in the 
number of correlators to be computed in order to achieve a prescribed minimal- match can be obtained using cardinal 
interpolation, thanks to the quasi-band-limited property of the match function [p9l , ]3(| . These results were obtained 
under the simplest assumption of a square-cell template-lattice. 

Two obvious questions are now in order: i) whether/to what extent cardinal interpolation is still effective when 
using the optimum triangular-tiling discussed in the previous sections |j3l| , and ii) whether/to what extent cardinal 
interpolation is still useful at relatively low T values, as needed in the early step(s) of hierarchical searches. Numerical 
simulations have been run to clarify both issues. 

In Fig. 8 the template density reduction p3| obtained by cardinal interpolation among optimally-placed templates 
is displayed as a function of the prescribed minimal-match T. In the inset of Fig. 8, the (boosted) minimal-match 
values obtained after cardinal interpolation among the optimally-placed templates corresponding to a number of 
representative values of T are also listed. 

At r = 0.97 a template density reduction sa 2.79 is obtained. The apparent discrepancy between this value and the 
one reported in ]30[ | is readily explained. The template density required to achieve a prescribed minimal-match when 
using cardinal interpolation is essentially independent from the shape of the lattice fundamental-domain, depending 
only on its measure (area), according to the well-known Nyquist-rate condition |33| ]. 

Thus, the template density reduction w 4 in ]3Q] , JlTj when using the simplest square-cell lattice is seen to result 
from two factors: a factor sa 2.79, expressing the template density reduction due to cardinal interpolation when using 
the optimum triangular-tiling, times a factor sa 1.43 expressing the template density reduction implied in using the 
optimum triangular-tiling instead of the simplest square one. 



V. CONCLUSIONS 



We presented a simple and systematic procedure for constructing the sparsest lattice of templates subject to a given 
minimal- match constraint in the range between r ~ 1 down to T sa 0.7 (and below) in the (spin- free) Tanaka-Tagoshi 
parameter-space of post-lPN gravitational wave chirps. 

We also showed that cardinal interpolation can be effective for correlator-bank economization both in the late and 
early stage(s) of multi-step hierarchical searches. 
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CAPTIONS TO THE FIGURES 



Fig 1 - Match contour-lines in Tanaka-Tagoshi coordinates. Curvature p 1 of 7r vs. polar angle <j> (LIGO-I, 2.5PN). 



Fig. 2 - Triangular, square and hexagonal tilings and pertinent Voronoi sets. 



Fig. 3 - The ratios r p , p = 3,4,6 relevant to Eq. as functions of V (LIGO-I, 2.5PN). 



Fig. 4 - Inscribed triangle with largest area does not necessarily cope with minimal-match condition for r < T* 
(LIGO-I, 2.5PN, L = 0.9). 



Fig. 6 - Optimum triangular-tilings and template-lattices for several values of T (LIGO-I, 2.5PN). 



Fig. 7 - Lattice-cell area corresponding to optimum triangular-tiling and number of templates needed to cover the 
range 0.2M Q < mi < m 2 < 1OM vs. minimal-match (LIGO-I, 2.5PN). 



Fig. 8 - Cardinal-interpolation gain (template density reduction factor and minimal-match boost) for optimum 
triangular-tiling (LIGO-I, 2.5PN). 



Fig. 5 - The butterfly-shaped set [PGQ 1 U QGP' n S r (G) (LIGO-I, 2.5PN, L = 0.875). 




Fig. 1 - Match contour-lines in Tanaka-Tagoshi coordinates. 
Curvature p~' of y r vs. polar angle </> (Ligo-I, 2.5 PN). 




Fig. 2 - Triangular, square and hexagonal tilings and pertinent Voronoi sets. 
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3 - The ratios r p , p = 3, 4, 6, relevant to Eq. (III.2), as functions of T 

(LIGO-I, 2.5PN). 




Fig. 4 - Inscribed triangle with largest area does not necessarily cope 
with minimal-match condition for T < r*(LIGO-I, 2.5 PN, T=0.9). 




Fig. 5 - The butterfly- shaped set (PGQ' vjQGP' )nS r (G) (LIGO-I, 2.5PN, r=.875). 





Fig. 6 - Optimum triangular-tilings and template-lattices 
for several values of T (LIGO-I, 2.5PN). 




Fig. 7 - Lattice-cell area corresponding to optimum triangular-tiling 

and number of templates needed to cover the range 
0.2M© < m, < m 2 < 10M© vs. minimal match (LIGO-I, 2.5PN). 
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Fig. 8 - Cardinal interpolation gain (template density reduction factor 
and minimal-match boost) for optimum triangular tiling (LIGO-I, 2.5PN). 



